Ultrasensitive Proximity Josephson Sensor with Kinetic Inductance Read-Out 
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We propose a mesoscopic kinetic-inductance radiation detector based on a long superconductor- 
normal metal-superconductor Josephson junction. The operation of this proximity Josephson sensor 
(PJS) relies on large kinetic inductance variations under irradiation due to the exponential temper- 
ature dependence of the critical current. Coupled with a dc SQUID readout, the PJS is able to 
provide a signal to noise (S /N) ratio up to ~ 10 3 in the THz regime if operated as calorimeter, while 
electrical noise equivalent power (NEP) as low as ~ 7 x 10 -20 W/\/Hz at 200 mK can be achieved 
in the bolometer operation. The high performance together with the ease of fabrication make this 
structure attractive as an ultrasensitive cryogenic detector of THz electromagnetic radiation. 

PACS numbers: 85.25.Oj, 74.45. +c,73.23.-b,73.50.Lw 



Superconducting single-photon detectors [J, 0, Q of- 
fer high infrared detection efficiency, high-speed timing 
resolution and few-nanosecond reset times. They have 
been applied in several fields including spectroscopy of 
ultrafast quantum phenomena (JL optical communica- 
tions quantum cryptography [y], and fast digital cir- 
cuit testing Q. On the other hand, a wide potential for 
superconducting nanoscale detectors used as advanced 
bolometers is also expected in several astrophysical space 
applications, where bolometers are promising candidates 
to meet future needs of cooled telescopes. The inter- 
est lies in the negligible Johnson noise they show with 
a NEP as low as 1CP 18 W/VHz- Hot-electron resistive 
microbolometers and kinetic inductance superconducting 
detectors (KIDs) represent high performance devices able 
to reach NEP< ICT 19 W/VHz at T > 1 K [8]. KIDs @ 
offer about the same NEP and response time as resis- 
tive bolometers and hot electron detectors, and they can 
operate at temperatures much below the critical temper- 
ature where the generation-recombination noise is small 
thanks to the reduced number of quasiparticles. 

Here we propose a KID based on a long 
superconductor-normal metal-superconductor Josephson 
junction. It exploits large kinetic inductance variations 
under irradiation thanks to the exponential temperature 
dependence of the supercurrent, and yields a high S/N 
ratio (~ 10 3 around 40 THz) and a low NEP (~ 7 x ICT 20 
W/VHz at 200 mK). The ease of implementation com- 
bined with large array scalability make this structure 
promising as a sub-Kelvin ultrasensitive detector of far- 
and mid-infrared electromagnetic radiation. 

The structure we envision is sketched in Fig. 1(a) and 
consists of a diffusive normal metal (N) wire of length I 
coupled to two superconducting leads (S) through trans- 
parent contacts, thus realizing a SNS Josephson weak 
link. An antenna eventually couples the incident radi- 
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FIG. 1: (color online) (a) Scheme of the PJS. Incident electro- 
magnetic radiation elevates the electron temperature (T e ) in 
the N wire, thus strongly suppressing the Josephson current. 
This leads to a large enhancement of the junction kinetic in- 
ductance, (b) Supercurrent 7 C vs electron temperature T e of 
a long SNS Josephson junction calculated at <f> = 7r/2 from 
Eq. Q (full line), and for A/E Th = 230 (dashed line). The 
inset shows a sketch of time evolution of T e after absorption 
of radiation (see text). 



ation to the N wire. We assume that the Josephson 
junction is long, i.e., A ^> HD/l 2 = E^h, where A 
is the S gap, D is the diffusion coefficient of N, and 
£/rh is the Thouless energy. The radiation coupled to 
the junction heats the electrons in N to temperature 
T c . For i?Th "C fcflT c <C A, the Josephson current is 
Ij = 7 c sin((/>), where 6 is the phase difference across su- 
perconductors, and [10( 
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is the junction critical current. Hence, in this limit, I c 
depends exponentially on the electron temperature and is 
independent of the phonon temperature Tbath- In Eq. |T]), 
i?N = pi I A is the normal-state resistance of the junction, 
p = (v-pe 2 D)~ 1 is the wire resistivity, A its cross section, 
and v-p is the density of states at the Fermi level in N. For 
our simulation we choose a 10-nm-thick silver (Ag) wire 
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with I — 1/im, width of 100 nm (volume 17 = 10 21 m 3 ), 
z/p = 1.0 x 10 47 J^m" 3 , and D = 0.01 mV 1 . With 
the aforementioned parameters i?N ~ 38 f2, and E^h — 
6.6 /xeV. By choosing, for instance, Nb as S electrodes 
(A = 1.52 meV) we get A/E Th ~ 230, thus providing the 
frame of the long junction limit. The critical current I c vs 
T is shown in Fig. 1(b) for (f> = vr/2 and A/E Th = 230 
(dashed line). In our case, I c saturates around 1.7 /xA 
at T c ~ 50 mK, and is suppressed by a factor of ~ 20 
at 1 K due to the exponential dependence on T c . For 
a comparison, the approximated result from Eq. (TTJ) is 
also shown (full line), and will be used in the following. 
Such suppression of I c produces a large enhancement of 
the junction kinetic inductance (L k ), defined as L k = 
h/(2el c ). As we shall show, measuring L k variations with 
a suitable readout scheme allows to accurately detect the 
radiation absorbed by the SNS junction. 

In order to understand the operation principle of the 
PJS as a calorimeter (i.e., in pulsed excitation operation) 
as well as a bolometer (i.e., in continuous excitation op- 
eration) it is useful to consider the inset of Fig. 1(b), 
which shows a sketch of time evolution of T e in N after 
the arrival of a photon. We assume that T c is elevated 
with respect to Tb a th, depending on the energy of the im- 
pinging photon and uniformly along N, over a time scale 
set by the diffusion time td = l 2 /D [see the red line in 
the inset of Fig. 1(b)]. With our parameters td = 10~ 10 
s. Then, after the absorption of a photon, T c relaxes to- 
ward Tbath over a time scale set by the electron-phonon 
interaction time (r c _ ph ), given by T c _ ph = l/(aT b 3 ath ) PJJ, 
where a « 0.34£/(fcg^F), and £ is the electron-phonon 
coupling constant. By setting S = 5 x 10 8 Wm -3 K -5 , as 
appropriate for Ag [11| . r c _ p h ~ lx 10~ 4 . . . 1 x 10 -7 s in 
the 0.1 ... 1.0 K temperature range, so that r c _ p h ^> td 
[see the blue line in the inset of Fig. 1(b)]. In the fol- 
lowing we analyze both the pulsed and the continuous 
excitation mode of the PJS. 

In the pulsed mode, after the arrival of a photon of 
frequency v at time t = 0, the electron temperature 
in N can be determined by solving the heat equation 
C c (dT c /dt) = P opt [IH, where C e = (7r^ F fc|T c )/3 is 
the electron heat capacity, and P op t = (2nhu/fl)S(t) is 
the optical input power per volume per incident pho- 
ton. In writing the heat equation we neglected the spa- 
tial dependence of T G in N, as well as the interaction with 
the lattice phonons, the latter occurring on a time scale 
Te-ph 3* tq. From the solution of the heat equation we 
get T c (v) — \/T 2 ath + I2hv / {itSlvFk^) , which shows that 
small N volumes are required to achieve large enhance- 
ment of T c . This condition can be easily met in metal- 
lic SNS junctions, where N island with volumes below 
10~ 21 m 3 can be routinely fabricated with the present 
technology. The relative variation of the kinetic induc- 
tance, i.e., the quantity SL^/L^. = [L^{v) — Lk(0)]/Lk(0) 
is displayed in Fig. 2(a) as a function of v at different 
Tbath- In the present structure, SL^/L^ of about 14% 
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FIG. 2: (color online) (a) 5L k /L|J vs v at different T ba th- (b) 
SL^/L^ vs P op t calculated at different Tb at h- 



for 1 THz photon, and around 163% for 10 THz photon 
can be achieved at 200 mK. At higher bath temperatures 
SL k /L^ is reduced, and obtains values of about 3% at 1 
THz and around 35% at 10 THz at T bath = 1 K. Such 
kinetic inductance variations allow for a very large signal 
to noise ratio for single-photon detection. 

The PJS operation in continuous excitation can be 
described by considering those mechanisms which drive 
power into the N electrons. At low temperature (typi- 
cally below 1 K), the main contribution in metals is re- 
lated to electron-phonon heat flux which can be mod- 
eled by Qc-ph = £ft(T 5 - T b 5 ath ) [llj. The steady- 
state T e under irradiation with a continuous power P opt 
thus follows directly from the solution of the energy bal- 
ance equation P opt + Qc- P h = 0, which gives T e (P opt ) = 



\J (Popt/Sfi) + Tb ath . This expression shows that both 
reduced 11 and small S are required to maximize T e 
enhancement upon power irradiation. The impact of 
continuous power excitation on the junction kinetic in- 
ductance is shown in Fig. [2jb) which shows SL^/L^. = 
[ik(Popt) - T k (0)]/L k (0) versus P opt at several T ba th. 
Notably, 5L k /L° as large as ~ 130% for P opt = 10 fW 
and ~ 2000% for P opt = 1 pW at T ba th = 0.2 K can 
be achieved. At higher bath temperatures SL^/L^. gets 
reduced, reaching values ~ 1% for P opt = 10 fW and 
~ 100% for P opt = 1 pW at 1 K. 

We now turn on discussing the PJS performance by 
considering a superconducting quantum interference de- 
vice (SQUID) readout, as shown in Fig. 3(a) [12] ■ A 
constant bias current p, divides into two parts, i.e., one 
flowing through the SNS junction [13], and the other (pj 
through a load inductor (L) coupled to a dc SQUID. 
Upon irradiation, an enhancement of L k results in a vari- 
ation of 7l , thus producing a magnetic field which is de- 
tected by the SQUID. The magnetic flux generated by 
the incident radiation is given by $ = MJl, where M 
is the mutual inductance between the SQUID and the 
SNS junction loop. In the linearized regime, i.e., by as- 
suming L/l <C 'I'o where $o is the flux quantum, we get 
h w / b #o/(*o+iic), and dl h /dl c w LI h <S> Q / (<S> +LI C ) 2 . 
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FIG. 3: (color online) (a) Scheme of the KID with a dc SQUID 
readout, (b) S/N ratio vs v calculated at different Tbath- (c) 
Temperature dependencies of NEP due to thermal fluctuation 
noise (TFN) and to the SQUID readout, (d) Resolving power 
vs v calculated at different Tb at h- 



In the pulsed detection mode, the signal to noise ratio 
(S/N) can be readily expressed as 



S_ 

N 



(d<P/dT c )5T c 



M(dI L /dI c )(dI c /dT c )6T c 



(2) 



where (5$ n is the flux sensitivity of the dc SQUID, and 
lu its bandwidth. The S/N ratio versus v is shown in 
Fig. 3(b) at different T ba th- Here we set L = 100 nH, 
M = 10 nH, lu = 1 MHz, 5$ n = 10- 7 $ /v / Hz Q, 
and lb = 0.8/ c (^). Notably, very high S/N ratios can 
be achieved with the PJS in the 100 GHz - 100 THz 
frequency range. The S/N ratio is maximized around 40 
THz where it obtains values ~ 1.2 x 10 3 at Tbath = 0.2 K. 
In the bolometer operation, on the other hand, an impor- 
tant figure of merit is the NEP, which is due to several un- 
correlated noise sources. In our case, the dominant con- 
tribution is due to thermal fluctuation noise- limited NEP 
(NEPtfn), given by NEP TFN = V^W+lLJ 
while the contribution due to Johnson noise is ab- 
sent, thanks to the operation of the junction in the dissi- 
pationless regime. The contribution of the SQUID read- 
out to NEP (NEPsquid) can be determined by setting 
S/N = 1, w = 1 Hz, and solving Eq. © for P opt . Figure 
3(c) shows the NEPtfn (dashed line) and NEPsquid 
(full line) vs Tbath- NEPsquid is significantly smaller 
than NEPtfn, and the latter can be as low as ~ 7x 10~ 20 
W/a/Hz at 0.2 K. Further reduction of NEPtfn is possi- 
ble by lowering f2 as well as by exploiting materials with 
lower E. Above, we have discussed the electrical NEP 
- the optical NEP is of the same order of magnitude. 
This is because the resistance of the device can be easily 
matched to common broadband self-similar lithographic 
antennas. The PJS resolving power (2nhv/ ' AE) vs fre- 
quency, where AE ss 2\/21n2NEPTFN(^)- v /''c-ph is the en- 



ergy resolution of full width at half maximum [ll| , is dis- 
played in Fig. 3(d) for different Tbath- In particular, the 
figure shows that resolving power values between ~ 1.2 
and ~ 2.3 can be achieved in the 5 ... 70 THz frequency 
range for Tbath ^ 400 mK, thus making the PJS suitable 
for far- and mid-infrared single-photon detection. 

The mechanism for the supercurrent in SNS junctions 
is due to the proximity effect, giving rise in the N lo- 
cal density of states to an energy minigap of size E g = 
c(</>)£rh with c(0) « 3.1 [TH, which we have ignored in 
the expressions for the heat capacity and electron-phonon 
coupling. Due to the minigap, both of these quantities 
are reduced inside the N wire, further improving the de- 
vice resolution. Moreover, the density of states is not 
divergent at the minigap edge, contrary to bulk super- 
conductors, and thereby the generation- recombination 
noise is reduced. For the chosen parameters the mini- 
gap Eg ~ h ■ 5 GHz, and thereby radiation from low 
frequencies v < E g /h, which typically is not part of the 
measured signal, should not couple in the device. 

In summary, we have analyzed a proximity Joscphson 
sensor (PJS) based on a long SNS junction in the kinetic 
inductance mode. Signal-to-noise ratio as high as ~ 10 3 
and noise equivalent power below 10 -19 W/vHz at 0.2 
K have been found to be achievable. Together with the 
available resolving power, the PJS is a promising candi- 
date for single-photon detection in the THz regime. 
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